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ABSTRACT 



M n thnH int g nHnd fnr A method providing r eal-time determination of a thermal 
profile of the drilling fluid in a well from three measuring points available in the field, 
i.e. the injection temperature, the outlet temperature and the bottomhole temperaturejs 
disclosed . The form of the profile between these three points is defined by a type curve 
representative of the thermal profiles in a well under drilling, estimated from physical 
considerations on thermal transfers in the well. 




FTFT .D OF TT - TR INVENTION Field of the Invention 

rOOOH The present invention relates to a method of determining the thermal 
profile of a drilling fluid in a well. 

[00021 During drilling, which is the mud injected into the drill string of the well 
and which flows flowing b ack through the corresponding annulus i undergoes great 
temperature variations. The fluid can encounter temperatures that can range from 
2°C for deep offshore wells to more than 180°C for very hot wells. Many mud 
properties, such as rheology or density, depend on the temperature. Calculation of 
the pressure losses during drilling can therefore be improved if an estimation of the 
temperature profile in the well is known. It is therefore important to be able to 
predict the temperature profile in the flowing mud from well data and mud 
characteristics. 

r00031 Measurement of the thermal profile of the fluid in a well under drilling 
would require complete instrumentation of the well, ire^ that is installation of evenly 
spaced out detectors in the drill string and in the annulus, allowing temperature 
measurement at various depths. However, installing such a measuring system 
entails too many constraints ; only localized measurements picked up by devices 
mounted in the drill string allow to know knowing certain temperature points on the 
path of the drilling fluid. 



2 



BACKGROUND OF THE INVENTION 

r00041 In the face of this lack of data, analytic models based on heat transfer 
equations have been developed to evaluate the thermal profiles of the fluid along the 
well under drilling. Some of these analytic models are implemented in softwares 
and allow to provide providing an estimation of thermal profiles from a certain 
numbor of data mnm nr latin which may be difficult to obtain. Thus, knowing the 
characteristics of the site and of the drilling equipment, by giving a value of the 
temperature of the fluid at the well inlet, these softwares can predict the temperature 
profile of the drilling fluid. 

f00051 However, a comparison between the results obtained with analytic 
methods and the measurements obtained in the field shows that there can be great 
differences. Furthermore, the complexity of softwares using numerical calculation 
methods makes real-time implementation thereof difficult. 

r00061 O n the other hand, a study of the bibliography on thermal models shows 
a similarity in the form of temperature profiles in most cases, which turns on three 
points : inlet temperature, outlet temperature and bottomhole temperature. 

r00071 The aim of this study is thus to propos e provide a method allowing real- 
time determination of a thermal profile in the mud from three measuring points 
available in the field, irer- that is the injection temperature, the outlet temperature and 
the bottomhole temperature measured by a detector mounted on the string. The form 
of the profile between these three points is represented by a type curve 
representative of the thermal profiles in a well under drilling, estimated from 
physical considerations on thermal transfers in the well. 
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SUMMARY OF THE INVENTION 

r00081 The method of determining the thermal profile of a drilling fluid 
circulating in a well under drilling according to the invention is defined by the 
successive stages as follows : 

a) determining a general expression 01 for the thermal profile of the fluid within 
the drill string in the well and a general expression 02 for a thermal profile of the 
fluid in the corresponding annulus, using the heat propagation equation that takes 
into account a thermal profile of the medium surrounding the well, 

b) measuring the temperature of the fluid at the well inlet, Tl, at the well bottom, 
T2, and at the well outlet, T3, 

c) laying down that modifying the expressions 01 and 02 to meet the temperature 
boundary conditions Tl, T2 and T3, 

d) drawing determining the thermal profile of the drilling fluid as a function of the 
depth. 

r00091 In order to obtain, in real time, a temperature profile with the method 
presented above, stages b), c) and d) can be repeated. 

F00101 According to the method of the invention, in stage a), general 
expressions 01 and 02 can comprise unknown constants, and in stage c), it can be 
laid down determined that expressions 01 and 02 meet the temperature boundary 
conditions Tl, T2 and T3 by determining said-the unknown constants. 

r00111 In order to determine a general expression 01 for the thermal profile of 
the fluid within the drill string in the well and a general expression 02 for a thermal 
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profile of the fluid in the corresponding annulus, it is possible, according to the 
method of the invention, to use in stage a) the heat propagation equation that takes 
into account at least the thermal equation of the medium surrounding the well, the 
flow rate of the fluid and the balance of the thermal exchanges undergone by the 
fluid, said-the thermal exchanges comprising at least the exchanges between the 
ascending and descending drilling fluid, and/or to use the equation of heat 
propagation in a homogeneous medium on a cylinder of infinite height centered on 
the well, said-the cylinder comprising the drill string that guides the descending 
fluid and the annulus around said drill string, which guides the ascending fluid. 

1*00121 According to the method of the invention, general expressions 01 and 92 
obtained in stage a) can be split up into several independent equations and, in stage 
c), it can furthermore be laid down determined that the profiles and the derivatives 
of the thermal profiles of the fluid within the drill string and in the corresponding 
annulus are continuous. 

r00131 The method according to the invention can notably be used to calculate 
the pressure drops of the drilling fluid circulating in a well under drilling, or in 
another application, to determine the zones of hydrate formation in the fluid during 
drilling. 

r00141 In relation to the methods for determining the thermal profile of a 
drilling fluid in a well according to the prior art, the present invention notably 
affords the following advantages : 

- the temperature profile obtained is more accurate because it is determined from 
three drilling fluid temperature measurement points while keeping an analytic 
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expression of the thermal profile between the measuring points which is physically 
justified, 

- since the temperature measurements being are performed all the time, the 
method allows to obtain obtaining the temperature profile in real time and to 
observe the evolution thereof with time. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[00151 Other features and advantages of the present invention will be clear from 
reading the description hereafter of non-limitative examples, with reference to the 
accompanying drawings wherein : 

. Figure Fig. 1 diagrammatically shows the architecture of a well under drilling, 

. Figur e s Figs. 2, 3 and 4 show the form of the temperature profile of the drilling 
fluid in a vertical onshore well, 

- Figure Fig. 5 shows the form of the temperature profile of the drilling fluid in a 
vertical offshore well, 

. Figure Fig. 6 shows the form of the temperature profile of the drilling fluid in a 
deviated offshore well, 

. Figur e Fig. 7 shows the evolution as a function of time of the temperature profile 
of the drilling fluid in a vertical offshore well. 



DETAILED DESCRIPTION 
DETAILED DESCRIPTION OF THE PREFERRED EBMODIMENTS 

r00161 It is possible to give an analytic expression for the thermal profile in the 
well and the annulus by using quite simple heat exchange considerations,4re r that is 
the heat propagation equation. 

1001-71 This model is based on the establishment of the heat balances in the 
well. According to a first approach, only the steady states are considered (the 
drilling mud flow is assumed to be stabilized for some time so that the temperatures 
no longer evolve). Certain hypotheses are necessary for calculation : the heat 
exchanges are measured in a plane perpendicular to the laminar flow of the mud, the 
various constants are assumed to be independent of the temperature, and finally the 
influence of the temperature of the medium surrounding the well shows on an 
apriori selected useful diameter Rf. 

r00181 It is then sufficient to use the heat propagation equation in a 
homogeneous medium on a cylinder of infinite height centered on the well shown in 
Figure 1 . In each well section, the heat loss equality is written by considering two 
temperature functions : G^z) within the drill string and 0 2 (z) in the annulus. 

TO0 191 Let 

0 f be the temperature of the formation, 

X f the thermal conductivity of the medium surrounding the well, 
X* the thermal conductivity of the tubing (metal), 
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Cp the heat-capacity rate of the drilling fluid, 
Rl the inside radius of the drill string, 
R2 the outside radius of the drill string, 
Rt the radius of the annulus, 

Rf the effective radius (for heat supply) around the well, 
D the flow rate of the drilling fluid, 
p the density of the drilling fluid. 

[00201 The heat balances per unit of depth are as follows : 

- Heat supplied by the medium surrounding the well to the fluid in the annulus 

in 1 f 
Qi = —? — f -x(02 



In 



V R f J 



Heat carried from the fluid in the annulus to the fluid within the drill string 



In 

- Heat accumulated by the fluid in the drill string and in the annulus 

Q, = -D.p.CpAG, 
Q a = D.p.C p A9 2 . 

r00211 The heat balances lead to the following system : 

Q. = Q 2 

Qa = Q. + Q 2 , 
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ire/That is 



d0 2 
dz 



dG, 
dz 



2TU f — 2m « (0,-* 2 ) 



DpC\n 



DpCAn 



2ru„ 



r00221 T hese equations are solved by diagonalization and matrix inversion , and 
they lead to the following results : 



0,(2) = -K x Be rXz - K 2 Be r2z +G f - — 

B 



6 2 {z) = + r,)e r -K 2 (B + r 2 )e r2z +G f 



with : 



A = 



2ru ( 



DpC p In 



r. = 



A + Ja 2 +4AB 



B = 



DpC p In 



\ R \ J 



A-«J A 2 + 4AB 



9 f = a.z+0o being the thermal equation of the medium surrounding the well and a 
the thermal gradient. 

K, and K 2 are the integration constants depending on the boundary conditions. 

r00231 It is thus possible, by means of some simplifying hypotheses, to obtain 
an analytic expression for the temperature profile of the drilling fluid in a well. If all 
the parameters are known, by giving the inlet temperature and by writing that 
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temperatures 9, and 0 2 are equal at the well bottom, the profile is entirely 
determined. The most commonly used softwares use this type of predictive 
procedure. However, a study of the results of the models compared to field data 
shows how difficult it is to use these estimations in a predictive way. 

r00241 In the present invention, the system is based on the knowledge of three 
measuring points in the field : inlet temperature, outlet temperature and bottomhole 
temperature. In order to estimate the thermal profile in the well from the three 
measurements consisting of the surface injection and outlet temperatures and the 
bottomhole temperature (inside or outside the drill string), the method according to 
the invention consists in connecting connects the three measuring points by a 
general expression representative of the evolution of a thermal profile in a wellbore, 
as obtained according to the method described above. 

F00251 Wo therefore Therefore take-the equations obtained by means of heat 
exchange calculations as follows are used : 

^(z) = -K x Be rX 2 -K 2 Be r2z +0 f 

0 2 (z) = -K X (B + r x )e rX z - K 2 (B + r 2 )e r2 z + 9 f 

r00261 According to the invention, these curve forms are adjusted to the three 
measuring points of the drilling fluid temperature at the inlet, Tl, at the well 
bottom, T2, and at the well outlet, T3. In order to use these three measuring points 
as boundary conditions, we choose to decoupl e the two equations (in the drill string 
and in the annulus) are decoupled by using different integration constants while 
keeping the general expression. Wo obtain two Two g eneral expressions of the 
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temperature profile are obtained in the drill string, 01, and in the annulus, 02, which 
have a physical significance but which comprise two degrees of freedom. Thus, 
expressions 01 and 02 can be adjusted by fixing said-the_degrees of freedom in 
order to meet the temperature conditions Tl, T2 and T3. Wo thorofor o decide that 
Therefore the equations in the pipes and in the annulus have the form as follows : 

0 t (z) = -K x Be rXz - K 2 Be r2z +e f ~ 

9 2 {z) = -K 3 (B + r x )e rXz -K 4 (B + r 2 )e r2z + 6 f . 

r00271 Wo thus havo four Four_integration constants K,, K 2 , K 3 and K4 are 
obtained rather than two, which requires four boundary conditions to determine the 
temperature profile. These four boundary conditions are : measurements of the well 
inlet, bottomhole and outlet temperature, and an equality condition at the bottom 
between the temperature in the drill string, 01, and the temperature in the annulus, 
02. The profile can be adjusted at any time to pass through the measuring points : 
wo thus have therefore a real-time estimation of the thermal profile is obtained . 
Programming by means of a spreadsheet type software allows to readily obtain 
ready obtaining of the real-time evolution of the profile. 

[00281 Figur e s Fjgs^_2, 3 and 4 respectively show the temperature profile of the 
drilling fluid in a vertical onshore well at a flow rate of 500 1/min, 1000 1/min and 
2000 1/min. The analytic expression determined allows simple calculation of the 
temperature T in Celsius degrees of the fluid in the drill string (curve 01) and in the 
annulus (curve 02) as a function of the depth P in meter. The analytic expression 
depends on several parameters that can be fixed from the start. Wo uso by By 
default typical values of these parameters are used . To determine the temperature 
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profile of Figures Figs. 2, 3 and 4, the geothermal gradient a is assumed to be 
constant to correspond to the onshore situation of the well. The temperature profile 
is entirely determined by measuring the temperature, 20°C at the inlet, 35°C at the 
bottom and 24°C at the outlet of the well. 

' f00291 The case of the vertical offshore well can be tackl e d addressed by 
considering that the geothermal profile of the medium surrounding the well is 
divided in two domains : let 9m be the thermal profile of the sea and 9s the thermal 
profile of the ground. Thermal gradient a is assumed to be constant in each domain, 
but discontinuous exist from one domain to the other. Let am be the thermal 
gradient of the sea and as the thermal gradient of the ground. Wo then consid e r two 
Two series of equations (one for each domain) are considered for each general 
expression in the pipes and in the annulus. We thus obtain four Four decoupled 
equations are obtained which represent the thermal profile of the drilling fluid in the 
well. Equation 91 l(z) corresponds to the temperature profile in the drill string in the 
sea, 912(z) corresponds to the temperature profile in the drill string in the ground, 
921(z) corresponds to the temperature profile in the annulus in the ground and 
922(z) corresponds to the temperature profile in the annulus in the sea, 911 being 
independent of 912 and 921 being independent of 922 : 



0 n (z) = -K x Be r ' 2 -K 2 Be' 22 +e m - — 

On 00 = ~K 3 Be^ - K 4 Be' 2 > + 0 M ~ 

0 2l (z) = -K 5 (B + r x )e rx '* - K 6 (B + r 2 )e r2z + 0 5 
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r00301 This brings the number of integration constants to eight (K, to Kg). The 
boundary conditions are then : measurements of the inlet, outlet and bottomhole 
temperature, equality condition at the bottom between the pipe temperature and the 
annulus temperature, to which wo add the continuity of the thermal profiles is added 
in the drill string and in the annulus at the junction of the two domains and the 
continuity of the derivative of the thermal profiles in the drill string and in the 
annulus at the junction of the two domains. Similarly, it is then possible to obtain in 
real time a physically realistic thermal profile passing through the measuring points. 
Figttre- Fig. 5 shows the thermal temperature profile of a drilling fluid in an offshore 
well from the four equations 911, 912, 921 and 922. The fluid circulates at 
509 1/min and the temperatures measured are 29°C at the inlet, 15°C at the outlet 
and 30°C at the well bottom. The thermal gradients are selected constant in each 
domain crossed by the well. 

r00311 Deviated wells represent the majority of the current wellbores. The 
physical problem is not fundamentally different and it can be handled in the same 
way as offshore wellbores : the well just has to be divided into two domains, each 
domain being characterized by a different thermal gradient corresponding to the 
medium surrounding the well. In the case of a deviated well, the depth corresponds 
to the distance covered along the well trajectory. General expressions 91 and 92 
representative of the thermal profile are each split up into two independent 
equations. The vertical part is characterized by the thermal gradient a of the 
medium surrounding the well, the deviated part is characterized by an equation of 
the thermal profile of the medium surrounding the well 9 d = cc.sin(<|>).z+9o, <|» being 
the angle of inclination. The same boundary conditions (measurements of the 
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temperature at the inlet, at the outlet and at the bottom of the well, equality at the 
bottom between the pipe temperature and the annulus temperature, and continuity of 
the thermal profiles and of the derivative of the thermal profiles in the drill string 
and in the annulus at the junction of the two domains) then allow to solv e solving 
the equations and to obtain obtaining the expression of the temperature profile in 
the pipes and in the annulus. 

r00321 It is possible to combine the procedure applied for a vertical offshore 
well and the procedure applied for a deviated onshore well in order to determine the 
temperature profile in an offshore well when the direction of the hole is deviated in 
the ground. The domain is divided into three different domains : let 0m be the 
thermal profile of the vertical domain in the sea, 9s the thermal profile of the 
vertical domain in the ground and 9d the thermal profile of the deviated domain in 
the ground. Fkwe- Fig. 6 shows the thermal profile in a deviated offshore well. The 
fluid circulates at 500 1/min and the temperatures measured are 20°C at the inlet, 
23°C at the bottom and 15°C at the outlet of the well. 

r00331 According to the same method as that used for the vertical offshore well 
or for the deviated onshore well, it is possible to determine the thermal profile of a 
vertical onshore well whose formation thermal gradient changes as a function of the 
depth. The well is divided into domains characterized by a thermal equation of the 
medium surrounding the well. General expressions 91 and 92 representative of the 
thermal profile are then each split up into as many independent equations as there 
are different domains. The same boundary conditions (measurements of the 
temperature at the inlet, at the outlet and at the bottom of the well, equality at the 
bottom between the pipe temperature and the annulus temperature, and continuity of 
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the thermal profiles and of the derivative of the thermal profiles in the drill string 
and in the annulus at the junction of the two domains) then allow to solv e solving 
the equations and to obtain obtaining the expression of the temperature profile in 
the pipes and in the annulus. 

r00341 B y repeating the calculation allowing to obtain obtaining the expression 
of the temperature profile of the drilling fluid upon each new temperature 
measurement, w o obtain a representation of the temperature profile evolving with 
time is obtained . Figure- Fig. 7 shows the evolution of the temperature profile of the 
drilling fluid in an offshore well in the course of time. The graph in the upper part 
of gjgttfe- Fi g. 7 shows the evolution as a function of time t in seconds of the flow 
rate parameter D in 1/min of the drilling fluid, and of the temperature parameter T in 
°C of the drilling fluid at the inlet, Tl, at the bottom, T2, and at the outlet, T3, of the 
well. The three graphs in the lower part show the temperature profile at three 
different times and allow to observe the evolution thereof. 

r00351 K nowledge of the thermal profile of the drilling fluid at any time allows 
real-time calculation of the pressure drops in the well by taking into account the 
thermal effects. This gives a better estimation of the bottomhole pressures and of 
the injection pressures for complex wells. 

r00361 Another use of real-time determination of the thermal profile of the 
drilling fluid is hydrates formation prevention. Hydrates form under low 
temperature and high pressure conditions, conditions which are met notably in deep 
offshore wells at the ground/sea interface. Knowledge of the temperature profile 
allows tn Hntmniino determination of t he zones where the temperature of the drilling 
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fluid is below the minimum value from which hydrates form, then to react 
accordingly, for example by raising the flow rate or by heating the fluid in order to 
prevent this formation of hydrates. 



